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A  SEA  ICE  TERRAIN  MODEL  AND  ITS  APPLICATION 
TO  SURFACE  VEHICLE  TRAFFICABILITY 

by 

W.D.  Hibler  III  and  S.P.  Ackley 


Utrodacttoa 

It  is  well  known  that  the  arctic  pack  ice  presents  a  formidable  surface  for  any  vehicle  traveling 
at  high  speed.  The  primary  obstacles  to  the  movement  of  amphibious  surface  vehicles  such  as  a 
surface  effect  vehicle  are  pressure  ridges,  which  reach  heights  as  great  as  20  ft,  and  often  are 
several  kilometers  long.  A  typical  example  of  a  pressure  ridge  is  shown  in  Figure  1,  which  gives 
a  ground  view  of  a  “first-year"  pressure  ridge  7  to  8  ft  high. 

The  statistical  occurrence  and  height  distribution  of  pressure  ridges  have  been  tho  subjects 
of  several  recent  studies  (Hibler  et  al.  1972,  Mock  et  al.  1972,  Weeks  et  al.  1971).  In  addition, 
considerable  progress  has  been  made  in  estimating  the  ridging  characteristics  of  the  Arctic  Basin, 
using  laser  profilometry  data  (Hibler  et  al.  in  prep).  For  classifying  pressure  ridges,  theoretical 
ridge  height  and  spacing  probability  density  functions  have  been  derived  which  have  been  found 
to  be  in  good  agreement  with  extensive  submarine  sonar  (Hibler  et  al.  1972)  and  laser  profilometer 
(Hibler  et  al.  in  prep)  data.  Using  these  theoretical  ridge  height  and  spacing  probability  density 
functions,  it  has  been  found  that  by  considering  ridges  to  be  composed  of  short  “ridge  links’*  the 


figure  1.  Ground  view  of  "first  year ”  pressure  ridge  in  the  Beaufort  Sea,  1971. 
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ridge  characteristics  of  a  given  region  may  be  well  defined  by  two  parameters:  the  mean  ridge 
height  h  and  the  ridge  density  RD,  which  is  defined  as  the  total  length  of  ridging  per  unit  area. 
Furthermore,  it  has  been  shown  (Mock  et  al.  1972)  that  for  randomly  oriented  ridgeB  fiD  is  given 
by  (n/2)fi  where  p  is  the  number  of  ridges  per  unit  length  intersected  by  a  straight-line  path.  Con¬ 
sequently,  by  using  digital  filtering  techniques,  F  and  RD  may  be  obtained  remotely  from  laser 
profile  data  (Hibler  1972).  The  final  information  needed  to  characterize  a  ridging  system  in  three 
dimensions  is  the  distribution  of  the  lateral  extents  of  ridges.  We  have  recently  obtained  estimates 
of  this  parameter;  the  results  are  reported  in  the  Appendix. 

Using  these  statistical  distribution  models,  a  realistic  three-dimensional  ridge  model  for  sea 
ice  terrain  was  constructed  in  this  study.  Using  this  model,  and  assuming  a  vehicle  height  clear¬ 
ance  ability,  a  trafficability  analysis  was  performed  by  calculating  the  ratio  of  the  total  distance 
traveled  Dr  to  the  component  of  distance  traveled  along  the  desired  travel  direction  °SL-  The  re¬ 
sulting  trafficability  ratios,  denoted  by  D T/DSL  as  well  as  other  parameters,  give  an  estimate  of 
transit  times  and  fuel  consumption  as  a  function  of  ridging  parameters  and  provide  c  basis  for  more 
detailed  mobility  analysis,  using  information  on  maneuverability,  obstacle  avoidance  systems,  and 
fuel  consumption  rates  as  a  function  of  speed.  Good  agreement  was  found  between  trafficability 
ratios  calculated  from  the  model  and  those  determined  from  simulated  routes  through  sea  ice  areas 
obtained  from  aerial  photo  mosaics. 

For  application  of  this  sea  ice  terrain  model  and  trafficability  analysis,  the  regional  variation 
in  ridging  intensity  in  the  Arctic  Basin,  based  upon  analysis  of  laser  profilometer  data,  is  discussed 
in  this  report.  In  particular,  the  number  of  ridges  per  kilometer  above  a  given  height  for  different 
regions  and  different  heights  is  discussed.  This  is  a  key  parameter  in  the  trafficability  model  and 
allows  estimation  of  the  ease  or  difficulty  of  motion  of  vehicles  having  a  given  clearance  ability  in 
different  regions  of  the  Arctic. 

Finally,  certain  other  aspects  and  ramifications  of  the  pack  ice  terrain  relating  to  amphibious 
vehicle  mobility  are  discussed,  with  emphasis  on:  1)  the  nature  of  the  shear  zone  formed  in  the 
offshore  region,  which  forms  a  rugged  but  relatively  thin  barrier  to  travel;  and  2)  the  presence  of 
long  linear  leads  which  occur  during  dilation  of  the  pack  ice  in  strong  weather  systems  and  which 
form  effective  high-speed,  flat,  parallel  "highways"  for  amphibious  vehicle  travel. 


Previous  Work 

Most  trafficability  and  mobility  studies  over  sea  ice  terrain  to  date  have  treated  the  terrain  in 
a  simplified  manner  by  considering  it  to  be  an  isotropic  two-dimensional  grid  of  point  obstacles. 
Typically,  the  differential  height  from  one  point  to  another  is  calculated,  using  a  multivariate  dis¬ 
tribution  function  based  upon  spectral  densities  of  the  surface  roughness.  Examples  of  such  studies 
are  a  detailed  treatment  by  Smith  and  Nakano  (1973),  and  less  detailed  discussions  ,by  Grumman 
Aerospace  Corp.  (1972)  and  LeSchack  and  Long  (1972). 

Although  useful  for  Monte  Carlo  calculations,  such  models  appear  to  bear  little  resemblance  to 
an  actual  sea  ice  terrain.  The  reasons  for  this  are  twofold.  The  first  and  basic  reason  is  that,  al¬ 
though  sea  ice  terrain  may  be  isotropic  on  a  large  scale  (say  30  km),  it  is  certainly  not  isotropic 
locally  because  ridges  are  generally  linear  features  sometimes  extending  several  kilometers. 
Consequently,  once  an  impassable  ridge  is  encountered,  the  vehicle  must  go  around  it  to  proceed. 
Simply  stated,  the  ridge  may  not  be  treated  as  a  point  obstacle. 

The  second  reason  is  that  an  isotropic,  multivariate,  normal  distribution  function  does  not 
provide  an  adequate  description  of  sea  ice  terrain.  Although  such  a  description  is  useful  for 
isotropic  Gaussian  stationary  terrain,  it  is  not  useful  for  sea  ice,  which  appears  to  be  very  aniso¬ 
tropic  on  the  local  level,  with  a  distribution  of  ice  heights  much  closer  to  lognormal  (Hibler  aod 
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Mock  1071)  than  Gaussian.  Certainly  the  local  anisotropy,  as  discussed  above,  is  a  problem  be* 
cause,  although  an  autocovariance  along  a  given  straight*line  path  across  the  terrain  may  indicate 
no  correlation  between  points  100  m  apart,  points  along  the  crest  of  a  ridge  persist  in  height  for 
hundreds  of  meters. 

Apart  from  isotropy  problems,  the  probability  of  obstacle  encounters  calculated  by  such  func¬ 
tions  appears  physically  unrealistic.  Smith  and  Nakano  (1973),  fa  example,  estimate  horn  a 
typical  spectral  density  that  a  2*ft*high  obstacle  occurs  with  a  0.0005  probability  over  a  40*ft 
interval,  a  about  one  2-tMUgh  obstacle  every  8  miles.  Since  laser  profile  data  indicate  that  a  2-ft- 
high  obstacle  typically  occurs  at  spacings  of  300  ft,  the  Smith  and  Nakano  predictions  are  certainly 
unrepresentative.  Therefore,  as  a  rule,  although  obstacle  heights  might  be  calculated  from  spectral 
densities  by  making  certain  assumptions  (Cartwright  1962),  such  uses  of  sea  ice  spectral  densities 
appear  unjustified  at  present. 

A  more  realistic  way  of  producing  a  terrain  model,  which  has  been  discussed  by  Aerojet- 
General  Cap.  (1972),  is  to  construct  random  ridge  heights  and  spacings  generated  accading  to 
height  and  spacing  distribution  models  [developed  by  Hibler  et  al.  (1972)1  along  a  given  path,  and 
then  to  calculate  the  average  detour  lengths  using  an  average  ridge  extent.  Such  model  calcula* 
tions,  although  of  considerable  usefulness  for  low  ridge  densities,  are  not  valid  for  high  ridge 
densities  because  the  model  does  not  take  into  account  ridge  intersections.  Fa  high  ridge  densi¬ 
ties,  there  are  a  large  number  of  such  ridge  intersections  and  cul-de-sacs  which  make  the  detours 
very  long  and  which  must  therefore  be  taken  into  account. 

The  model  developed  in  this  study  not  only  takes  into  account  intersecting  ridges  but  alsq 
recognizes  the  peculiarities  of  the  sampling  problem;  fa  example,  given  ridges  of  equal  length,  a 
vehicle  selectively  encounters  the  ridges  nearly  prependicular  to  its  path  more  often  than  it  does 
those  approximately  parallel  to  its  path.  Also,  given  ridges  of  similar  orientation,  but  different 
lengths,  a  vehicle  encounters  the  longer  ridges  mae  frequently  than  it  does  the  shorter  ones.  In  a 
complete  trafficability  calculation,  these  factors  must  be  included  and,  as  is  shown  later,  they  do 
have  important  ramifications. 


Sea  Ice  Tetraia  Model 

The  basic  defect  with  most  sea  ice  trafficability  studies  to  date  has  been  the  unrealistic  sea 
ice  terrain  models  used.  Consequently,  the  emphaslB  in  this  study  has  been  to  develop  a  realistic 
three-dimensional  sea  ice  obstacle  model  which,  although  making  trafficability  calculations  more 
difficult,  is  a  reasonable  representation  of  the  actual  terrain.  To  create  such  a  model,  we  assume 
the  following  characteristics  based  upon  ridge  statistical  studies. 

1.  The  distribution  of  pressure  ridge  heights  above  some  arbitrary  cutoff  height  hQ  is  given  by 
the  probability  P  density  function  (Hibler  et  al.  1972): 


P(h)dh 


2/A\*  exp  (-to*)dh 
\nJ  afc  (\/XhQ) 


(1) 


where 

P(A)dh  =  the  fraction  of  ridges  with  heights  between  ft  and  ft  +  dft 
erfc  =  the  complementary  error  function 

A  =  the  shape  parameter  which  is  determined  uniquely  by  5,  the  mean  ridge  height. 
[Shape  parameters  fa  different  regions  of  the  Arctic  Basin  are  given  by  Hibler  et  al.  (in  prep).i 
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2.  The  probability  at  occurraoce  at  ridges  along  a  linear  path  is  given  by  a  Poisson  distribu¬ 
tion  (Hibler  et  al.  1972).  Consequently,  the  spacing  distribution  is  a  negative  exponential  and 
both  it  and  the  occurrence  distribution  are  well  described  by  fi,  the  average  number  at  ridges  per 
kilometer. 

3.  The  orientation  of  ridges  is  assumed  to  be  random  (Mock  et  al.  1972);  that  is,  for  a  given 
area  there  are  approximately  equal  lengths  of  ridges  oriented  at  any  angle.  As  a  corollary  to  this, 
it  can  be  shown  that  the  ridge  density  RD  (total  length  of  ridges  per  square  kilometer)  is  related  to 
H  along  a  linear  path  by  RD  =  (n/2)fi. 

4.  The  linear  extent  of  a  ridge  before  an  appreciable  hole  occurs  in  it  is  described  by  some 
distribution  function  PL(L)dL,  with  first  and  second  moments  <L>  and  <LS>.  This  function  must 
be  estimated  from  aerial  photo  mosaics;  and  as  a  first  approximation,  we  aesume  that  Pl(L)  is 
approximately  independent  of  ridge  height  for  ridges  higher  than  5  ft  and  is  negative  exponential  in 
form.  An  estimation  of  Pt(L)  from  pressure  ridge  data  is  given  in  the  Appendix. 

5.  All  the  distribution  functions  are  assumed  to  be  independent;  for  example,  if  a  high  ridge 
is  crossed  there  is  no  preferred  probability  that  the  next  ridge  crossed  will  also  be  high  (or  low). 
This  ha?  been  tested  to  some  extent  (Hibler  and  Mock  1971). 

Based  on  these  characteristics,  a  typical  sea  ice  obstacle  model  could  be  created  in  the 
following  conceptual  manner.  First  a  “box'*  of  random  ridges  could  be  created  using  the  height 
distribution  P(h)dh  and  the  length  distribution  Pl(L)dL.  Then,  taking  a  two-dimensional  region, 
a  number  of  these  random  (random  both  in  x  and  y  coordinates  and  in  orientation)  ridges  would  be 
placed  on  the  surface  until  the  total  length  of  ridges  equaled  the  desired  ridge  density  RD. 

Certain  qualifications  to  the  above  model  should  be  noted.  One  qualification  is  that  a  profile 
along  a  ridge  rapidly  varies  in  height.  Consequently,  although  the  use  of  the  lateral  extent  distri¬ 
bution  is  reasonably  valid  for  larger  vehicles,  small  vehicles  such  as  snowmobiles  may  be  able  to 
slip  through  a  ridge  at  many  locations  because  of  small  holes  in  the  ridge.  This  is  discussed  in 
somewhat  more  detail  by  Hibler  and  Ackley  (1973). 

Another  point  is  the  variation  of  PL(L)dL  for  different  ridge  heights.  Since  such  variations 
may  occur,  trafficability  analyses  should  generally  provide  for  determining  the  sensitivity  to 
changes  in  PL(L)dL.  The  assumption  of  random  orientation  is  alBO  one  that,  although  reasonable 
over  large  areas  (tens  of  kilometers),  breaks  down  over  small  areas  (several  kilometers)  where 
ridges  often  seem  to  occur  in  a  linear  pattern  (Hibler  1971,  Mock  et  al  1972).  The  height  and 
spacing  distributions,  however,  do  appear  to  be  valid  for  all  ridging  systems. 


Trafficability  Model 

To  carry  out  a  simple  trafficability  calculation,  we  will  consider  a  completely  maneuverable 
vehicle  that  can  clear  any  ridge  up  to  a  height  h.  We  will  then  calculate  the  ratio  of  the  total 
distance  traveled  to  the  component  of  distance  in  a  desired  direction.  This  ratio  £>t/DSl  turns  out 
to  be  a  unique  function  of  the  number  of  ridges  per  unit  length  higher  than  h,  denoted  by  and  of 
the  first  and  second  moments  (<L>  and  <L2>)  of  the  ridge  extent  distribution  function.  Conse¬ 
quently,  using  statistics  compiled  from  laser  profile  data,  the  Dj/Dsl  ratios  may  be  estimated  for 
different  regions  of  the  Arctic. 

First,  an  analytic  calculation  with  ridge  intersection  neglected  is  performed.  This  analytic 
calculation  is  similar  to  that  discussed  by  Bekker  (1969)  with  extensions  to  include  obstacles  of 
different  lengths.  This  calculation  is  useful  for  illustrating  how  the  sampling  biasing  affects  the 
results;  i.e.,  long  ridges  perpendicular  to  the  desired  motion  are  encountered  more  frequently  than 
short  ones.  Then  the  basic  calculation  for  DT/D SL  s  performed,  using  a  Monte  Carlo  method  which 
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takes  into  account  cul-de-sacs  and  ridge  intersections.  The  point  at  which  the  results  obtained 
with  the  Monte  Carlo  method  approach  those  obtained  with  the  analytic  model  is  also  important 
since  it  indicates  at  what  point  ridge  intersection  problems  become  negligible.  In  general,  as 
might  be  expected,  the  results  indicate  that  DT/DgL  very  rapidly  increases  as  ridge  intersections 
begin  to  be  encountered. 

Analytic  calculation  of  *V°«L 

Consider  an  area  of  dimensions  W  by  W  with  N  ridges  higher  than  h  dispersed  randomly  through¬ 
out  in  both  orientation  and  direction.  The  lengths  of  the  ridges  are  such  that  there  is  a  probability 
PL(L)dL  of  a  ridge  having  a  length  between  L  and  L  +  dL.  We  consider  that  a  desired  vehicle  path 
begins  at  the  bottom  of  the  area  and  proceeds  north.  We  neglect  lateral  offset  because  of  the  random 
orientation  of  ridges.  The  probability  that  any  given  ridge  intersects  the  vehicle’s  path  is  given  by: 

LP,(L)dL 

— — I***!"  (2) 

where  0  (— w/2  <>  9  z  n/2)  Is  the  angle  of  the  ridge  relative  to  north.  For  each  ridge  of  length  L 
and  orientation  0,  the  vehicle,  on  the  average,  travels  an  additional  distance  L/2  (1  -  cos  0)  going 
around  the  ridge.  Summing  over  all  ridges  in  the  area  and  over  all  angles,  we  have  the  extra  dis¬ 
tance  AO: 


AD-4"  /  //2  (1  -  cos  0)  |  sin0 1  id 

0  -n/2 

N  o 

=  — —  vLz>  /  (sin0  -  sin0  co*-.  0)d0 


N<LZ> 

2nW  ' 

However,  the  total  number  of  ridges  intersecting  the  vehicle  path  is: 


N  °°  .  ~  «r/8 

Ph"  /  ! sin 0 1 d0 

n  0  -n/2 

2N<L> 

=  ~~WH  ' 


Therefore,  the  total  distance  traveled  over  the  straight-line  distance  is 


W  +  AD 
V 


N<LZ> 
- +  1 

2vW* 


1  <L8> 

4  <L> 


Ph  +  1 


(3) 


(4) 


(5) 


where  ^  is  the  number  of  ridges  per  unit  length  with  height  greater  than  h  encountered  along  a 
straight-line  path. 
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The  most  striking  aspect  at  the  analytic  result  is  that  the  trafficability  ratio  Dr/DgL  depends 
an  the  distribution  of  ridge  lengths  only  through  the  ratio  of  the  second  moment  to  the  first  moment. 
Also,  the  factor  <LZ>/<L>  only  appears  linearly,  so  that  small  changes  in  Pl(L)  are  not 
extrenely  critical  in  determining  mobility;  this  is  certainly  an  advantage  since  Pl(L)  is  a  difficult 
distribution  function  to  estimate  exactly  for  sea  ice. 

Monte  Carlo  calcnlation 

Clearly  the  primary  defect  in  the  analytic  approach  is  the  neglect  of  ridge  intersections.  In 
order  to  detfmine  the  effect  of  such  intersections  on  trafficability,  as  well  as  to  verify  the  analytic 
model,  a  Monte  Carlo  calculation  was  performed.  The  calculation  proceeded  by  the  following  steps. 

1.  A  set  of  random  ridges  that  would  be  encountered  along  a  straight-line  path  was  constructed. 
The  number  of  random  ridges  having  a  given  orientation  0  and  length  L  were  weighted  by  a  factor 

L  coed  where  0  is  the  angle  between  the  ridge  and  the  straight-line  path.  This  factor  takes  into 
account  the  incieased  probability  of  encountering  longer  ridges  nearly  perpendicular  to  the  path  of 
motion.  The  distribution  of  ridge  lengths  was  taken  to  be  a  negative  exponential  and  was  con¬ 
structed  so  that  <LZ>/<L>  =  1.6  km. 

2.  Using  small  increments  AX,  we  proceeded  along  a  desired  straight  line,  searching  randomly 
for  a  ridge  in  each  increment  according  to  an  input  ridge  frequency  ^  (number  of  ridges  per  kilometer 
above  height  h). 

3.  If  no  ridge  was  encountered,  the  desired  distance  and  total  distance  traveled  "'ere  increased 
by  AX. 

4.  If  a  ridge  was  encountered,  the  vehicle  heading  was  changed  to  be  parallel  to  the  random 
ridge  of  length  L  at  direction  0  so  that  the  vehicle  still  had  a  forward  component  of  motion. 

5.  The  point  at  which  the  ridge  was  encountered  (and  thus  the  distance  to  go  around  the  ridge) 
was  randomly  selected. 

6.  We  then  proceeded  in  AX  increments  at  a  heading  0,  recording  after  each  increment  a  total 
distance  of  AX  and  a  desired  distance  of  (cos0)AX.  In  each  inciement,  a  random  check  for  other 
ridges  was  made. 

7.  If  the  end  of  the  ridge  was  reached  before  another  ridge  was  encountered,  the  heading  was 
changed  to  0=0  and  we  returned  to  step  2  and  proceeded  as  before. 

8.  If  another  ridge  at  direction  O'  was  encountered  while  <M0  (i.e.  while  we. were  still  going 
around  the  first  ridge),  we  declared  a  cul-de-sac  and  changed  the  heading  of  the  vehicle  to  direction 
O'  (or  -O')  so  that  the  vehicle  had  a  backward  component  of  motion.  We  then  proceeded  as  in  step 
five. 


9.  The  program  ended  once  the  total  distance  DT  reached  a  desired  value.  The  straight  line 
or  desired  distance  DSL  is  always  less  than  or  equal  to  f>T. 

10.  The  resulting  DT/DS,  ratio  and  number  of  cul-de-sacs  were  recorded  as  a  function  of  the 
dimensionless  parameter  ub<Lr>/<L>.  It  is  easy  to  see  that  DT/DgL  is  a  unique  function  of  this 
parameter  by  dimensions^  analysis. 


One  of  the  more  important  steps  in  the  Monte  Carlo  computation  is  step  8,  which  takes  into 
account  the  effects  of  intersecting  ridges.  The  procedure  used  does  not  account  exactly  for  such 
intersections,  but  simply  declares  that  the  vehicle  must  have  a  backward  component  of  motion  after 
encountering  a  ridge  intersection;  i.e.  it  must  proceed  along  a  ridge  so  that  it  has  a  negative  compo¬ 
nent  of  motion  in  the  desired  direction.  This  approximate  procedure  gives  a  slightly  high  estimate 
of  the  extra  distance  traveled  but  is  a  reasonable  first  approximation  commensurate  with  the  accuracy 
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Figure  2.  Traffic  ability  ratio  DT/DgL  as  a  function  of  ^  (<  L2  >)/(<  L  >). 

-  number  of  impassable  ridges/km,  <L>  -  first  moment  of  the  ridge 
length  distribution,  and  <L2>  -  second  moment  of  the  ridge  length  distri¬ 
bution. 

of  the  terrain  model.  Also,  In  practice,  since  the  extra  distance  traveled  rises  very  rapidly  as  more 
cul-de-sacs  are  encountered,  a  qualitative  estimate  of  the  extra  distance  is  adequate  for  most 
trafficability  purposes. 

Note  also  that,  similarly  to  the  analytic  result,  the  trafficability  ratio  Dr/DSL  depends  only 
on  <1  2>/<L>,  so  that  uncertainties  in  the  exact  moments  of  the  distribution  of  the  ridge  lengths 
may  not  be  critical. 

Res  alts  of  trafficability  compilations 

Random  simulations  were  made  using  the  Monte  Carlo  computer  program  described  above  with 
a  total  distance  for  each  simulation  of  100  km.  For  each  value  of  the  ridge  frequency  (and  hence 
each  value  of  the  parameter  ^  <L2>/<L>),  six  random  simulations  were  performed.  For  each 
simulation  the  trafficability  ratio  Dt/DSl  was  recorded,  and  average  results  of  the  six  independent 
simulations  were  computed  for  each  value  of  The  results,  as  a  function  of  <L2>/<L>,  are 
shown  in  Figure  2. 

As  can  be  seen  from  Figure  2,  the  effect  of  cul-de-sacs  forces  the  trafficability  ratio  up  rapidly. 
Since  typical  values  of  <L2'*/<L>  are  of  the  order  of  1.5  km,  the  results  indicate  that,  for  more 
than  2  impassable  ridges/km,  movement  becomes  rather '  jrtuour.  although  possible,  up  to  about  3 
ridges/km.  The  other  important  limiting  case  is  the  point  at  wnich  the  analytic  approximation 
approaches  the  Monte  Carlo  result.  From  f  igure  2,  this  would  appear  to  be  reasonably  valid  for 
** h  <L*  >/<L  >  <  l  or  for  the  presence  of  less  than  about  0.5  impassable  ridges/km. 

Experimeatal  trafficability  ratios  using  ridge  overlays 

In  order  to  test  the  trafficability  calculations  experimentally,  ridge  overlays  of  two  different 
regions  of  the  Arctic  were  prepared  using  aerial  photo  mosaics  with  typical  routing  paths  drawn 


Figure  4.  Rldgea  prepared  from  an  aerial  photo  mosaic  overlay  ot  sea  Ice.  Simulated  vehicle 
routtpga  are  Indicated  by  solid  lines  with  snows.  ,(Slte  oi  Match  1971  AIDJBX  camp.) 
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Figure  5.  Ridge  model  prepared  from  Figure  4,  with  the  addition  of  a  set  of  random  ridges  having 
the  same  length  distribution  as  those  shown  in  Figure  4.  Simulated  routing  paths  are  shown  by 

solid  lines  with  arrows. 

through  them  from  which  trafficabllity  ratios  were  calculated;  the  results  were  compared  with 
predicted  results  based  on  the  average  ridge  frequency  of  the  photo  mosaics.  The  overlays  were 
also  used  to  estimate  the  nature  of  the  distribution  oC  ridge  lengths  as  discussed  in  the  Appendix. 
The  overlays  were  prepared  by  using  photo  mosaics  from  which  dominant  ridges  were  outlined.  Al¬ 
though  there  was  considerable  approximation,  this  procedure  supplied  a  qualitative  measure  of 
ridges  above  6  ft  high,  especially  since  shadows  might  be  used  to  estimate  ridge  height. 

Figures  3  and  4  illustrate  two  ridge  overlays  with  typical  routing  paths  constructed.  The  ridge 
frequencies  were  0.31/km  for  Figure  3  and  0.65/km  for  Figure  4.  In  addition  to  these  two  photo 
mosaics,  a  photo  mosaic  having  a  higher  ridge  density  was  constructed  by  adding  a  set  of  random 
ridges  to  the  mosaic  in  Figure  4  having  the  sane  length  distribution  as  the  ridges  in  Figure  4. 

This  constructed  overlay,  with  a  ridge  frequency  of  1.06/km  is  shown  in  Figure  5. 

Besides  the  routings  drawn  on  the  photo  mosaics  in  Figure  5,  several  other  routings  were  also 
constructed  on  each  photo  mosaic  and  average  trafficabllity  ratios  were  computed.  A  comparison 
of  these  experimental  trafficabllity  ratios  with  predicted  ratios  based  on  the  Monte  Carlo  calcula¬ 
tion  and  our  estimated  values  of  ^  and  <Le>/<L>  is  given  in  Table  I.  The  predicted  and  observed 
results  are  in  reasonable  agreement,  with  the  observed  ratios  generally  being  slightly  smaller  then 
the  predicted  values.  Since,  as  discussed  earlier,  the  Monte  Carlo  calculation  tends  to  slightly 
overestimate  the  trafficabllity  ratio,  this  type  of  deviation  is  not  unexpected  and  might  be  particu¬ 
larly  pronounced  when  the  rlages  have  a  definite  lineation,  which  may  have  occurred  in  Figure  4. 


Tattle  1.  Comparison  of  predicted  and  observed  trafficabllity  ratios. 


Photo 

mosaic 

h 

<l8>/<l> 

<D^/DS^>  observed 

DTi'DSL  Pr0d<ct*d 
(Monte  Carlo) 

Fig.  3 

0.31/km 

1.88  km 

1.08  i  0.07 

1.10 

Fig.  4 

0.66/km 

1.64  km 

1. 16  i  0.09 

1.46 

Fig.  6 

1.06/km 

1.66  km 

1.71  ±  0.4 

1.80 
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Figure  6.  Vertical  aerial  view  of  sea  ice  showing  the  intersecting  nature  of  ridges.  Obtained 

the  Beaufort  Sea  by  NASA.  March  1971.. 
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These  ridge  overlays  also  give  a  perspective  oo  the  nature  of  terrain  for  different  ridge  fre¬ 
quencies.  Of  the  three  overlays.  Figure  5  with  a  ridge  frequency  of  about  1  ridge/km  shows  the 
most  dense  array  of  ridges  and,  as  is  apparv.r,  transit  across  such  a  region  would  not  be  rapid. 
However,  oo  the  scale  of  ridge  frequencies  encountered  in  the  Arctic,  as  discussed  in  the  next 
section,  there  are  often  more  than  5  ridges/km  above  4  ft  high  so  that  1  ridge/km  is  not  a  heavy 
ridging.  The  terrain  in  Figure  8,  on  the  other  hand,  with  a  ridge  frequency  of  about  0.8  ridges/km 
begins  to  approach  a  terrain  that  can  be  rapidly  moved  over,  although  even  in  this  case  some 
maneuvering  would  be  necessary.  Generally,  it  appears  that  for  rapid  transit  there  must  be  fewer 
than  0.5  impassable  ridge s/km. 

A  further  example  illustrating  the  nature  of  the  sea  ice  terrain  is  given  in  Figure  6,  which 
shows  a  vertical  serial  photo  of  ridging  taken  in  the  Beaufort  Sea  from  an  altitude  of  3500  ft.  This 
is  s  relatively  lightly  ridged  area  of  the  Arctic  Ocean,  as  discussed  in  the  next  section.  The 
photos  illustrate  again  the  intersecting  nature  of  the  ridging  and  reinforce  the  conclusion  that  ridge 
intersections  are  important  obstacles  for  vehicles  that  cannot  travel  over  most  of  the  ridges  en¬ 
countered. 


Regional  Variations  ia  Ridging  Intensity 

In  order  to  estimate  ease  of  travel  in  different  regions  of  the  Arctic  Basin,  using  the  above 
trafficability  analysis,  it  is  necessary  to  estimate  the  regional  variability  of  the  number  of 
ridges  above  height  A.  The  regional,  seasonal  and  annual  variations  of  ridging  intensity  have  been 
the  subject  of  considerable  study;  a  more  detailed  quantitative  discussion  of  such  variations  based 
upon  over  3000  km  of  laser  profile  data  taken  in  different  years  and  seasons  is  given  by  Hibler  et  al. 
(in  prep).  The  results  of  these  studies  of  variations  indicate  that  from  November  through  March  the 
overall  regional  ridging  variations  may  be  approximately  categorised.  Figure  7  shows  the  results  of 
such  a  categorisation  in  terms  of  the  number  of  ridges  above  given  heights  of  >4  ft,  >6  ft,  >8  ft  and 
>  10  ft.  Since  there  is  some  variation  in  mean  ridge  heights  from  year  to  year,  although  the  relative 
regional  variations  are  similar,  these  heights  may  have  about  a  1-ft  uncertainty.  For  a  year  of  heavy 
ridging,  for  example,  the  4-ft  category  refers  to  the  number  of  ridges  above  5  ft  and  for  a  year  of 
light  ridging  the  number  of  ridges  above  approximately  4  ft. 

As  can  be  seen  from  Figure  7,  the  western  portion  of  the  Arctic  Basin  is  divided  into  three 
ridging  sones  which  are  in  order  of  decreasing  ridge  heights:  1)  the  offshore  zone,  consisting  of 
the  region  off  the  Canadian  Archipelago  and  Greenland;  2)  the  central  arctic  zone;  and  8)  the 
Beaufort  Sea  zone.  Note  that  there  is  also  an  Alaskan  sooe,  which  includes  the  thin  region  off 
Alaska  and  the  Canadian  Northwest  Territory,  where  shearing  of  the  ice  pack  against  the  stationary 
shore  ice  produces  a  relatively  narrow  region  of  heavier  ridging  before  the  more  lightly  ridged 
Beaufort  Sea  zone  is  reached.  This  shear  sooe  tends  to  be  quite  variable  but  in  terms  of  ridge 
heights  seems  to  have  ridge  heights  similar  to  or  higher  than  those  of  the  offshore  sooe.  Other 
aspects  of  the  shear  sooe  are  discussed  in  the  next  section. 

There  is  some  difference  In  the  morphology  of  ridges  in  the  different  regions.  The  central 
arctic  sooe  has  more  multiyear  ridges  which  have  relatively  smooth,  rounded  surf  sees  than  the  other 
sones,  whereas  the  other  sooes  contain  many  first-year  ridges  which  an  rough,  often  consisting  of 
angular,  jumbled  blocks.  A  more  detailed  discussion  of  ridge  structure  s  is  given  by  Weeks  et  al. 
(1971).  To  convey  more  clearly  the  nature  of  different  ridges,  several  ground  views  of  various 
ridges  are  shown  in  Figures  8-10.  Figure  8  shows  a  multiyear  ridge  about  7  ft  high  which  is 
rather  smooth  and  would  be  relatively  easy  to  cross.  Figure  9  shows  a  small  first-year  ridge  about 
5  ft  high  and  Figure  10  shows  a  first-year  ridge  7-8  ft  high.  Clearly  the  rough  angular  nature  of  first- 
year  ridges  makes  them  more  difficult  to  clear  than  multiyear  ridges. 


Number  of  Ridge*  per  km  above  given  height 


Central 

Arctic 


Off  more 


4.0  to  7.0 
1.9  to4.0 
0.4  to  1.0 
0.2  to  0.5 
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Figure  11.  Trafficability  ratio  DT/Dgi  u  a  tuoctlon  ot  vehicle  clear¬ 
ance  ability  for  varlona  tonea  ot  the  Arctic  Baata.  The  aolld  curvet 
were  calculated  for  the  maximum  ridging  intensity  encountered  in  the 
various  regions.  The  dashed  curve  was  calculated  for  the  minimum 
ridging  intensity  encountered  in  the  Arctic  Basin  (Beaufort  Sea). 


It  is  clear  from  the  ridging  variations  shown  in  Figure  7  that  design  for  a  vehicle  operating  off 
Alaska  in  the  Beaufort  Sea,  and  perhaps  in  the  central  Arctic,  might  be  considerably  different  from 
that  for  a  vehicle  meant  to  operate  everywhere  in  the  Arctic.  A  vehicle  having  the  ability  to  clear 
5  to  6-ft  ridges,  for  example,  could  travel  reasonably  well  in  the  Beaufort  Sea,  but  would  have  con¬ 
siderably  decreased  mobility  in  other  zones,  including  the  shear  none  off  Alaska. 

The  effect  of  regional  ridging  variation  on  vehicle  trafficability  is  shown  more  clearly  in 
Figure  11,  which  shows  trafficability  ratios  as  a  function  of  vehicle  clearance  height  for  various 
ridging  intensities.  The  curves  are  obtained  by  using  typical  values  of  the  mean  ridge  height  and 
frequency  of  ridges  above  2  ft  hlgn  to  evaluate  the  height  distribution  function  given  in  eq  1. 

The  number  of  ridges  above  a  given  height  is  then  calculated  and  trafl.c ability  ratios  are  taken  from 
Figure  2  using  <LS>/<L>  =  1.5  km.  The  solid  curves  were  obtained  using  the  parameter  for  the 
mo6t  intense  ridging  found  from  laser  statistics  for  a  given  region.  The  dashed  line  was  calculatec 
using  parameters  obtained  for  the  least  intense  ridging  in  the  Beaufort  Sea.  The  curves  indicate 
the  sharp  increase  in  the  trafficability  ratio  for  a  given  region  as  the  height  clearance  is  decreased. 


Additional  Trafficability  Aspects  of  Boa  Ice 

There  are  several  further  aspects  of  the  arctic  pack  ice  that  relate  to  trafficability.  Two  such 
aspects  are  1)  the  shear  zone,  with  its  concomitant  rubble  and  hummock  fields,  where  the  pack  ice 
meets  the  stationary  shore-fast  ice;  and  2)  long  linear  leads  due  to  wind  driven  dilation  of  the 
ice  pack. 

Shear  sons  and  nibble  fields 

Often  in  the  shear  zone  constant  grinding  of  the  pack  ice  produces  a  dense  collection  of  rough 
ice  blocks  and  hummocks.  This  can  also  happen  in  localized  regions  elsewhere  in  the  Arctic.  An 
oblique  aerial  view  of  the  shear  zone  is  given  in  Figure  12.  In  extreme  cases  of  this  type  of  grinding 
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Figure  12.  Oblique  aerial  vie w  ot  the  shear  zone  oil  Barrow.  Alaska. 


Figure  13.  Closeup  view  ot  a  rubble  field  in  the  Beaufort  Sea. 


and  ridging,  as  discussed  by  Kovacs  (19711,  a  series  of  quite  high  (several  meters)  closely  spaced 
hummocks  or  ridges  forms  that  looks  like  a  farmer's  plowed  field.  In  less  intense  cases,  the  ice 
simply  forms  a  series  of  jumbled  blocks,  sometimes  with  large,  discontinuous,  vertical  faced  ridges, 
as  shown  in  Figure  13.  Figure  13  gives  a  closeup  of  a  rubble  field;  although  the  heights  of  the 
blocks  are  not  great  the  rough  surface  may  offer  considerable  hazard  to  certain  types  of  vehicles. 

A  descriptive  account  of  travel  through  this  zone  by  foot  is  given  by  Herbert  (1970). 

For  classification,  the  shear  zone  generally  is  expected  to  have  considerable  rubble  but  may 
not  always  have  extremely  high  ridges,  especially  not  at  all  locations.  Consequently,  as  a 
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trafficability  problem  tbe  shear  too©  can  be  considered  a  relatively  thin  but  formidable  barrier 
with  tbe  ridging  intensity  falling  off  farther  into  tbe  Arctic  Ocean  away  from  tbe  shear  sene.  For 
detailed  routing  studies,  the  shear  tooe  merits  considerably  more  study. 


As  opposed  to  most  aspects  of  sea  ice  terrain  discussed  above,  linear  lead  systems  tend  to 
aid  amphibious  travel,  at  least  in  ooe  direction,  rather  than  impede  it.  Studies  of  sea  ice  drift  and 
deformation  to  date  (Hibler  et  al.  1978,  Hibler  1978)  generally  indicate  that  during  winter  In  a  high 
pressure  region  tbe  ice  pack  tends  to  open  up.  A  typical  lead  system  present  under  such  a  dilation 
•  is  shown  in  Figure  14.  'Hie  larger  leads  shown  are  as  wide  as  500  m  and  extend  a  considerable 
distance.  The  availability  of  these  lead  systems  varies  with  the  weather  system,  and  Btudies  are 
presently  underway  to  predict  their  occurrence  from  meterological  data. 

Clearly,  if  ooe  were  Interested  in  traveling  parallel  to  a  lead  system,  much  more  rapid  vehicle 
progress  could  be  made  than  normal.  Also,  since  in  winter  such  systems  tend  to  form  in  high 
pressure  regions  with  scoompanylng  good  weather,  visibility  would  be  expected  to  be  better  than 
normal. 


Conclusions 

The  study  of  sea  ice  terrain  and  tbe  formulation  of  a  trafficability  model  indicate  that  inter* 
sections  of  pressure  ridges  make  travel  over  the  pack  ice  by  a  sizeable  (width  several  meters  or 
larger)  vehicle  difficult  when  the  vehicle  cannot  clear  most  pressure  ridges  encountered.  If  there 
are  more  than  abcut  0.5  Impassable  ridges/km,  rapid  travel  with  a  sizeable  vehicle  becomes 
difficult  and,  if  there  are  more  than  3  impassable  ridges/km,  any  sort  of  travel  with  a  sizeable 
vehicle  becomes  very  difficult. 

Using  these  numbers,  tbe  general  statistical  nature  of  arctic  pressure  ridges  indicates  that  a 
vehicle's  capability  of  clearing  5  to  6-ft  ridges  is  a  minimal  requirement  for  travel  in  the  Beaufort 
Sea  and  central  arctic  regions,  and  of  clearing  8-  to  9-ft  ridges  is  a  minimal  requirement  for  the 
offshore  zone  off  the  Canadian  Archipelago. 

In  addition  to  pressure  ridging  in  the  Arctic  Basin,  the  shear  zone  off  Alaska  and  the  Canadian 
Northwest  Territory,  with  its  rubble  and  hummock  fields  in  many  locations,  forms  a* thin  but  Very 
rugged  barrier,  although  ooe  could  probably  cross  this  barrier  in  a  vehicle  by  carefully  searching 
for  an  opening  in  the  barrier. 

One  possible  exception  to  the  trafficability  model  baaed  upon  ridges  may  occur  when  long 
linear  leads  form  during  dilation  of  tbe  pack  ice  in  high  pressure  regions  in  winter.  These  leads 
often  form  in  parallel  patterns  and  may  supply  rapid  transit  zones  in  ooe  direction  for  vehicles 
which  otherwise  would  be  farced  to  move  slowly  through  the  pressure  ridges. 
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APPENDIX  A.  DISTRIBUTION  OF  THE  LATERAL  EXTENT  OF  RIDGES 

Ridge  leaf  ;hs  are  a  somewhat  difficult  parameter  to  estimate  for  a  number  of  reasons.  One  of 
the  primary  reasons  is  that  the  nature  of  the  distribution  of  ridge  lengths  makes  it  necessary  to 
cample  a  relatively  .urge  area  (tens  to  hundreds  of  square  kilometers).  Since  ridges  are  relatively 
low  relief  features  (<  15  to  20  ft  high),  aerial  photography,  the  primary  source  of  data  for  this  type 
of  iniormation,  must  be  obtained  at  low  altitudes  (<10,000  ft)  in  order  to  adequately  resolve  these 
features.  However,  the  cameras  must  not  be  flown  so  low  that  the  area  covered  is  not  a  large 
enough  sample;  so  a  trade-off  between  obtaining  photography  that  adequately  resolves  ridges  and 
sampling  a  sufficiently  large  region  to  be  statistically  valid  must  be  achieved.  It  is  also  necessary 
to  use  photo  mosaics,  since  ridges  generally  extend  beyond  the  region  sampled  by  one  aerial  photo¬ 
graph.  This  constraint  also  limits  the  study  to  two  specific  areas  since  photo  mosaics  at  low  level 
on  sea  ice  are  not  generally  available. 

Once  adequate  data  are  available,  some  subjective  judgments  must  then  be  applied  to  obtain 
the  ridge  length  distributions.  The  first  concerns  the  definition  of  the  height  cutoff.  Low  ridges 
(< 3  ft)  are  generally  subdued  by  snowdrifting  and  other  relief  changes,  so  they  are  much  harder  to 
distinguish  than  high  ridges  (>  5  ft).  Also,  since  it  is  generally  agreed  that  vehicles  should  be 
able  to  traverse  'Tow”  objects,  and  the  purpose  of  this  study  was  to  obtain  a  length  distribution 
for  ridges  considered  dllticult  or  impossible  to  cross,  the  following  semiempirical  criteria  were 
used  to  define  ridges:  1)  Ridges  were  used  in  the  distribution  if  they  were  greater  than  -  5  feet 
high,  based  on  comparison  of  known  areas  from  ground  measurements  with  the  same  areas  portrayed 
on  the  aerial  photographs;  2)  Ridges  were  considered  to  end  when  the  height  level  dropped  into  a 
‘Tow”  category  (estimated  to  be  below  3  ft)  and  this  low  category  persisted  for  a  minimum  of  100  m 
(328  ft).  Each  continuous  segment  defined  in  this  manner  was  called  a  ridge  and  Us  length  was 
measured  at  the  appropriate  scale. 

Using  the  above  criteria,  ridge  overlays  of  two  photo  mosaics  were  prepared.  These  overlays 
are  shown  in  Figures  3  and  4  in  the  main  text.  The  distribution  of  the  ridge  lengths  and  other 
length  statistics  are  given  in  Figures  Al  and  A2.  As  can  be  seen,  the  values  of  <L2>/<L>  are 
of  the  order  of  1.2  to  1.5  km,  with  the  shape  of  the  distribution  tending  toward  a  negative  exponen¬ 
tial.  Fitting  a  negative  exponential  distribution  (of  the  form  N p  exp  (~px),  where  p  =  1/<L>  and  N 
is  the  number  of  ridges)  to  the  observed  distribution,  the  x2  totals  were  within  the  5%  confidence 
limits  in  both  cases.  A  negative  exponential  is  physically  reasonable  for  ridge  lengths  since  it 
can  be  derived,  for  example,  by  assuming  that  openings  in  ridges  occur  randomly  so  that  spacing 
between  holes  is  exponentially  distributed  (for  example,  see  Hibler  et  al.  1972).  We  note  in  passing 
that  for  a  negative  exponential  distribution  <L2  >/<L>  is  equal  to  2  <L> . 

For  generalisation  of  these  ridge  length  distributions  to  either  higher  (>.-  5  ft)  or  lower 
(<3  ft)  ridges,  examination  cf  available  aerial  photographs  and  field  observations  suggests  that 
the  mean  ridge  length  increases  slightly  with  increasing  height,  and  may,  for  example,  be  of  the 
order  of  2  km  fa  8  to  9-ft  ridges  and  1  km  fa  3  to  4-ft  ridges. 


Ptrcefrt  of  Ridgtt 
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Figure  Al.  Data  obtained  from  aerial  photo  mosaic  taken 

■  at  75°N,  145°W  in  the  Beaufort  Sea  [site  of  the  Marcb-April 

H  1972  AIDJEX  camp  (Fig.  3).]  Area  sampled  was  441  km* 

■  and  total  number  of  ridges  was  307.  Mean  ridge  length 

|ij  <L>  =  0.635  km.  Mean  aqua  re  length  <Le>  -  0.816  km2. 

Igl  Variance  ~  0.413  kr2.  <L2>  t <L>  -  1.28  km.  No.  of 

ridges/unit  length  (>  5  ft)  «  0.31/km. 
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Figure  A 2.  Data  obtained  from  aerial  photo  mosaic  taken 
at  74°N,  131°W  in  the  Beaufort  Sea  [site  of  the  March  1971 
AID J EX  camp  (Fig.  4}].  Area  sampled  was  332  km 2  and 
total  number  of  ridges  was  180.  Mean  ridge  length  <L>  « 
1.02  km.  Mean  square  length  <L*>  -  1.57  km2.  Variance 
»  0.536  km2.  < L2> /<L>  *  1.54  km.  No.  of  ridges/unit 
length  (>5  ft)  =  0.65/km.  Because  of  the  height  of  the 
photography,  ridge  lengths  shorter  than  0. 25  km  were  not 
included  in  the  distribution,  as  shown  in  the  figure. 
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